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We study the temperature dependence of single-particle and pair coherent condensate 
densities for a Bose liquid with a "depleted" single-particle Bose-Einstein condensate (BEC) 
at T 7^ 0. Our investigations were based on the field-theoretical Green's functions. On the 
basis of the use of empirical data about the speeds of the first and second sounds is described 
superfluid state at T 7^ 0. It is studied the structure of superfiuid state taking into account 
appearance with the different from zero temperatures of normal component and taking into 
account the branch of the second sound, whose speed approaches zero at T ^ Ta. The bare 
interaction between bosons was chosen in the form of a repulsive Azis potential, the Fourier 
component of which was an oscillatory and sign-varying function of the momentum. The 
obtained temperature dependence of single-particle and pair coherent condensate and total 
densities has a good agreement with experiment data. 
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1 Introduction 



Despite the big progress achieved in the theory of superfluidity since the pioneer works by Lan- 
dau PP, Bogolyubov 2|, Feynman [2| and others the task of constructing a microscopic 
theory of a superfluid (SF) state of a '^He Bose hquid cannot be considered complete. Analysis 
of experimental and theoretical activities concerning with unique superfluidity phenomena inves- 
tigation has one conclude that microscopic theory of superfluidity is still incomplete. There are 
some unsolved contradictions between theory and experiment as concerned with superfluid liquid 
hydrodynamics as with structure of quasiparticle spectrum. These contradictions reduce to the 
question of quantum-mechanical structure of superfluid "^He component below the A-point. This 
question is crucial in creation of consistent microscopic superfluid theory of Bose liquid both for 
the case of zero temperature and for the case of different from zero temperatures. 

In works one introduces at T = the microscopic model of the SF state of a Bose 

liquid with a single-particle Bose-Einstein condensate (BEC) suppressed because of interaction, 
based upon a renormalized field perturbation theory with combined variables |llj-jl4|. allows one 
to obtain a self-consistent "short" system of nonlinear integral equations for the self-energy parts 
T,ij{p, e), by means of truncating the infinite series in the small density of the BEC (ng/n <C 1). By 
the same token, one can work out a self-consistent microscopic theory of a superfluid Bose liquid 
and perform an ab initio calculation of the spectrum of elementary excitations E{p), starting from 
realistic models of pair interaction potential U(r). 

In this work we will study of the superfluid state structure of Bose-liquid *He under nonzero 
temperature (0 < T < Tx). The same task was consider in work 25 on the basic of Chester 
IBG wave function vIZ/bg by Jastrow factors F = Ufij. We use the microscopic model 
proposed in for the superfluid Bose liquid with a "depleted" BEC. The small parameter of 
these model is the ratio of the density of the BEC to the total density Bose liquid po/ ps 1 in 
contrast to the Bogolyubov theory |2] for a almost ideal Bose gas, where the small parameter is 
the ratio of the number of overcondensate excitations to the number of particles in the intense 
BEC (n - no) /no < 1. 

Such approach is based on the analysis of the numerous precise experimental data on neutron 
inelastic scattering [22], [2H|-|Sni and to results in quantum evaporation of ^He atoms ac- 
cording to which the maximal density po of the single-particle BEC in the ''He Bose liquid even at 
very low temperatures T <^T\ does not exceed 10% of the total density p of liquid *He, whereas 
the density of the SF component ps ^ p a.t T ^ Q ^U]. Such a low density of the BEC is im- 
plied by strong interaction between ^He atoms and is an indication of the fact that the quantum 
structure of the part of the SF condensate in He II carrying the "excess" density (pg — po) ^ po 
calls for a more thorough investigation. So, in this case the density of the superfluid component 
Ps is determined by the quantity of the renormalized anomalous self-energy parts Si2(0, 0), which 
is a superposition of the " depleted" single-particle BEC and the intense " Cooper" pair coherent 
condensate (PCC), with coincident phases (sings) of the corresponding order parameters. Such a 
PCC emerges due to an effective attraction between bosons in some regions of momentum space, 
which results from an oscillating sign-changing momentum dependence of the Fourier component 
V(p) of the interaction potentials U{r) with the inflection points in the radial dependence. 

The system of Dyson-Belyaev equations [S], which describe the superfluid (SF) state of a Bose 
liquid with strong interaction between bosons and a weak single-particle Bose-Einstein condensate 
(BEC) and which allows one to express the normal Gn and anomalous G12 renormalized single- 
particle boson Green functions in terms of the respective self-energy parts Sn and S12, has the 
next analytic form of |27j : 




(1) 
(2) 



where 



(3) 
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V{p, e) = V{p) [1 - y(p)n(p, e)]-' . (4) 

where V{p) is the Fourier component of the bare potential of the pair interaction of the bosons; 
V{j), e) is the renormahzed (screened) pair interaction between bosons because of many-particle 
collective effects; n(p, e) is the boson's polarization operator which describe many-particle corre- 
lation effects 

X Gii(k + p,e + tj) +G'i2(k,w)Gi2(k + p,e + tj)} ; 

where T(j),e]k^uj) is the vertex part (three-pole) describing many-particles correlations of the 
local field's type; A(p, e) = r(p, e, 0,0) — r(0,0,p, e) is the vertex part with zero values of input 
momentum and energy; and no is the number particles in BEC, ni = n — no is the number of 
overcondensate particles (ni S> no), which is determined by the condition of the total particle 
number conservation: 

f (Pk f duj ^ , , 

n = no + ni^no + i J j^-:^ J —Gii[k,uj). (6) 

Take into account the residues at the poles of single-particle Green functions Gij{p, e), neglecting 
the contributions of eventual poles of the functions r(p, e, k, lo) and V^(p, e), which do not coincide 
with the poles of Gij{p,e), the functions ^ij{p,e) on the "mass shell" e = E{p) assume the 
following form (at T — 0): 



vl/n(p,i?(p)) - 2 y j^r{p,E{p}-k,E{k)) 



X V{p-k,E{p)-Eik)) 

1 r d^k 



(27r)3 

A{KE{k)) 



(7) 



§i2(p,i?(p)) = -i J- 



(27r)3 

X Vip - k, Eip) - E{k))T{p, E{p); k, i?(fc)) (g) 
noA(k, E{k))Vik, E{k)) + *i2(k, E{k)) 



X 



where 



Aip) = noA(p, Eip))V{p, Eip)) + n.ViO) + *n(p) + ^ - (9) 
here /i is the chemical potential of the quasiparticles, which satisfies the Hugengoltz-Pines relation 

El 

M = Sn(0,0)-Ei2(0,0) . (10) 

The quasiparticle spectrum corresponding to the poles of Green's function is determined in general, 
with allowance for relations H1I2|) . by the next expression 



E{p) = ^A^p) - [noA(p, E{p))Vip, E{p)) + ^i2(p)] (11) 

2 The structure of the superfluid Bose hquid state at T 7^ 

The main purpose of our research is the investigation of superfluid helium condensate structure by 
means of finding temperature dependence of single-particle and pair coherent condensate densities 
under assumption that superfluid component is a superposition of a weak single-particle BEC and 
an intensive PCC. 
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Let's consider the superfluid bose liquid state at T ^ when there are superfluid component 
density ps (T) and normal component density p„ (T) simultaneously. As it was shown in [131 114j 
the expression for the renormalized Green's function Gij (p) is constructed with help of combined 
variables: 

^-(2;) = + (12) 

In the long- wavelength region |fc| < fco (where k is some characteristic momentum) this variables 
are just the hydrodynamics variables ^l(x) in the spirit of Landau quantum hydrodynamics 
and in the short-wavelength region |fc| > fco they coincide with the usual field operators 5's/i(a;): 



(13) 



^.h = ^ - V'l; Mr) = ^ E|k|<fco "ke^"^"" = Vii^e^-^- . 

Such an approach means that the separation of the Bose system into a macroscopic coherent 
condensate and a gas of supracondensate excitations is made not on the statistical level, like in 
the case of a weakly nonideal Bose gas but on the level of ab initio field operators, which are 
used to construct a microscopic theory of the Bose liquid. 

In the region of small p 7^ at T — s- Green's functions have the form 



1 Tl 

4-ff..(p)--2 



—9^-iP)-^'^Mp)--- (14) 



~ \ Ti ■ 

Guip) = risg^^ip) - —g^^{p) - Y'^vviP) ■ ■ ■ (15) 



where 



J\a\<an ^'^ 



'\q\<<ia 

p= {k,e),q= {q,u) (16) 

where g^ivip) are the "hydrodynamics" Green's functions, which are associated with the long- 
wavelength fluctuations of the phase and density of the condensate ( pLv = ipn). The expressions 
for gcptpip), gipn{p), gmrip) was calculated in |12) for T > 0. It contains sums of two pole terms, 
corresponding to first and second sound with velocities ci and C2 in the Bose liquid with the 
normal and superfluid components: 

/, X _ ja^u - df^^pn/p) b^^ pn/p 

p.,iy^ip,7T (17) 

where p = p„ + Ps is the total density of the liquid, and the coefficients a^i^^d^i, and are 
independent of T at low temperatures. This result establishes a unique correspondence between 
the microscopic field theory of superfluidity ISl [71 and the macroscopic two-fluid hydrodynamics 

It follows from Eqs. H14|) . IjlSfl and H17|l that the pole parts of the renormalized Green's 
functions Gij can be represented in the form 



We henceforth assume that expression H18|l is valid in the entire temperature interval T <T\. 
It follows from Eqs. lfTH|l that for T — > 0, where p„ 0, the leading contribution to the integral 
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over energy e in Q comes from the first-sound pole e = ci/c of the Green's functions. However, 
at higher temperatures T > IK, where p„ ~ ps, the main part, because of the strong inequahty 
ci ^ C2, is played by the low-energy pole e = C2k which corresponds to second sound. 

At finite temperatures (T ^ 0), after taking into account the contributions of the first and 
second poles of the Green's functions (|18|l . we obtain for the self-energy parts in F = A = 1.5 
approximation 



Tiij (fc, T) : 



2 



xcothf^)+i3,^^cothf^)l 
\2TJ ^ p C2q \2TJ\ 



1 
ciq 



(19) 



It should be emphasized that the long-wavelength approximation for the Green's functions 
((TH|l in this case are valid because of the divergence of the temperature factor coth(c2(;/2T) at 
g — > and the rather rapid decay of the interaction kernel ( q — *■ oo). Moreover, the system of 
equations (|19(l docs not need to be matched with the expression for the renormalized quasi-particle 
spectrum E(k), as it is ordinarily done in microscopic field theory for T — > 0. It's true because the 
substitution of the empirical spectra of the first and second sound (with the experimental values 
of the velocities ci and C2) into the expressions for the Green's functions Gij{p) corresponds to 
automatically taking all the necessary renormalizations into account. 

Using ()19|l . one can determine the superfluid order parameter for T 7^ 0: 



Si2(o,T) = *o(r)-f*,(r) 



Ps{T) 



(20) 



where 



*o(r)--- 



V{q) 



A12 - £'12 
ciq 



coth 



2TJ 



■ coth , 
C2q \2TJ 



(21) 



^^iT) = -'- [ fl,V{q) \^ coth f ^) - ^ coth f ^ 
^ ^ 2 J (27r)3 ^'[ciq \2T I C2q \2T 



(22) 



Since under T = the density of the superfluid component ps coincides with the total mass 
density of the Bose liquid p = mn and is proportionate to I]i2(0, 0) (it plays a part of superfluid 
order parameter) we can write the following expression 



Ps = Po + Ps ^ I3m 



Sl2(0) 

mv{o) 



/3m[no(l-7) + ^], 



(23) 



where 



P dk 



(27r)2A(0)y(0) 7 E{k) 



(A(fc)V^(fc) 



(24) 



(27r)2A(0)l/(0) 7 E{k) 



A(A:)l/(fc)*(fc), (25) 







where (3 is the constant to be identifled. Starting from the definition of the single-particle BEC 
density po as po = mno, we obtain /? = (1 — 7)"^. Thus the "Cooper" PCC density takes form: 

Ps = mni = , (26) 
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where po — mriQ is the density of the single-particle BEG and ps is the density of "Cooper" PCC. 

On the other hand, the temperature dependence of BEC density po{T) = mno{T) is due to 
entire particle number conservation condition ^ will be determined the following expression 



Po(T) 



1 

= 1 - - 
2 

X coth 



(2^ 

ciq\ 
2t) 



Pn{T) 1 



PniT) 



P 



1 

ciq 



P 



coth(f) 
C2q \21 / 



The density of pair coherent condensate is 

Ps{T) ^ vl/o(T) 
P V{Q)n 



n -1 



V{Q)n 



(27) 



(28) 



Thus the "Cooper" PCC density at T = in our model can be found in two independent ways. 
We can find it directly with help of expression 12()|) or we can subtract H27I) from (|28|l and then 
find result at T = 0. Thus we can verify the self-consistency of the proposed model via comparison 
the results of independent calculations of the PCC density of superfluid component for the zero 
temperature. 



3 The iterative scheme of calculation and results 

In order to calculate the the temperature dependence of single-particle and pair coherent conden- 
sate densities within the model of a Bose liquid with a suppressed BEC considered, at first, we 
have to find the Fourier component of the bare potential of the pair interaction of the bosons and 
then, to find the renormalized (screened) pair interaction between bosons, which is determined 
so we must to calculate the boson polarization operator 

The pair interaction between bosons was chosen in the form of a regularized repulsive potential 
in the the Aziz potential |82l : 



Aexp(-Q;r-/3r2) -exp [-(ro/r- 1)2] ^^^pC2fe+67- ^'^ ^ , r < tq 
Aexp(-ar - /Jr^) - Y.\=o C2k+(ir^'^''^^ , r > ro 



where ^ = 1.8443101 x lO"^ K, a = 10.43329537 A"!, ^ = 2.27965105 A"^^ cg = 1.36745214 Kx A , 
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C8 = 0.42123807 Kx A , cio 0.17473318 Kx A . Such potential remain finite at r = due to the 
nonanalytic exponential dependence on r, which suppresses any power divergence at r — s- 0. It has 
the Fourier component V{p), which is an oscillatory and sign- varying function of the momentum 
p as a result of the "excluded volume" effect and the quantum diffraction of the particles on one 
another. 

Then, in order to calculate the polarization operator, using Eqs. (U)) and ((SJ, a numerical 
calculation in the first approximation of the functions ^i{p) = ^ii(p, £'o(p)) and 'I'i(p) = 
^i2(p, i?o(p)), was conducted. For the zcroth approximation, the "screened" potential were 
taken: 

Voip) - , (30) 

pa - VoUoji (pa) 

at some constant negative value of Hq. Then, using the functions ^i{p) and 5*1 (p) obtained, 
the first approximation for the polarization operator Ili{p) was calculated, using Eqs. lfn|l - lfTC|l . 
(|18|l at r — 1. The limiting value ni(0) was compared with the exact thermodynamic value 
of the polarization operator of the "^He Bose liquid at p = and w = 0, which determines the 

71 

compressibility of the Bose system j^: 11(0, 0) = 5- . 

mcj 

The absolute value |n(0,0)| turned out to be almost 1.5 times greater than the calculated 
value |ni(0)|. This provides for an estimate of the vertex F at p = in the first approximation as 
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Fi = Ai ~ 1.5. The second approximation ^2{p) and ^'2(p) was obtained from Eqs. iQ, © with 
the constant value Fi = Ai and the first approximation for the renormahzed pseudopotential, 
Eq. gj): 

Such an iterative procedure was repeated four to six times and used to improve precision in 
the calculation of the polarization operator. At each stage, equations and were used 
to reproduce the quasiparticle spectrum E{p), and the rate of convergence of the iterations was 
watched, as well as the degree of proximity of E{p) to the empirical spectrum Ecxp{p)- The fitting 
parameter in these calculations was the amplitude Vq of the seed potential (|29|l at the value of 
a — 2.44 A, which is equal to twice the quantum radius of the ''He atom. The BEC density, in 
accordance with the experimental data [2HI~in3> was fixed at uq — 9%n = 1.95 • 10^^ cm~'^. 

For the numerical computation of temperature dependencies of single-particle and pair co- 
herent condensate density we have to take into account that the first (hydrodynamics) sound 
velocity is practically independent of T and in the given approximation can be determined as 

Ci = l/(0)n/m* ( here m* is the effective mass of quasiparticles) , whereas the velocity of 

second sound C2 is substantially T-dependent, varying from C2(0) = ci/\/3 at T = to a value 
C2(T) ~ 20m/s in the region T > IK 8 . At T ^ Tx the velocity C2 ^ 0. Thus as the A point 
is approached the main part, owing to the strong inequality ci 3> C2, begins to be played by the 
last terms in the integrands in (|21|l . H27|l . which are proportional to B12 and Bn and contain the 
temperature factor 

f{q,T) = ^—coth(^^)=^^, C2q<T (32) 
C2{T)q V 2fci3T J c4{T)q^ 

which diverges quadratically as g ^ 0. At the same time the width of the singular peak increases 
rapidly with increasing T and decreasing 02- (the momentum dependence of temperature factor, 
obtained for different temperatures, is depicted in fig. 1). As a result, with increasing T there is an 
increase of the contribution to the integral (|21|l from the repulsive part of the potential V{q) > 
in the long- wavelength region q < n/a and a decrease of the function '^o{T)- The function '^o{T) 
plays a part of the superfluid order parameter and is positive at low T < C2 g owing to the strong 
attraction V{q) < in the region ir/a < q < 27r/a. At a certain critical temperature T ^ Tc 
the function '^o{T) goes to zero and then becomes negative (for T > Tc), that corresponds to 
destruction of the superfluid state {ps = 0), i.e. Tc coincides with the A point. 

In a similar way with increasing T, there is increase of the negative contribution to the integral 
(|27|l and decrease of po{T) until at a certain point T = Tq the density of the BEC vanishes and 
formally becomes negative for T > Tq. 

The results of the numerical calculations of the temperature dependence of the densities of 
single-particle and paired parts of superfluid density in the range of temperatures from zero to the 
point of lambda-transition, are depicted on fig. 2. 



4 Conclusions 

Thus, on the basis of the use of empirical data about the speeds of the first and second sounds 
in this work is described superfluid state at T 7^ 0. It is studied the structure of superfluid state 
taking into account appearance with the different from zero temperatures of normal component and 
taking into account the branch of the second sound, whose speed approaches zero at T ^ Tx. We 
obtained the analytic expression 1)271 128|l for computation po and ps densities, and constructed and 
realized the numerical scheme for computation of the temperature dependence of single-particle 
and pair coherent condensate densities. The temperature dependence of the densities of single- 
particle condensate and total density of superfluid component are obtained on the basis of the 
microscopic model of superfluid state of Bose-liquid with the depressed single-particle BEC in the 
range from T = to temperatures close to the environment of point of A-transition. 
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Note also that the self-consistency of this model is corroborated by the fact that the theoretical 
value of total particle density calculated from Eq. @, nth = 2.1 • 10^^ cm~^, is quite close to the 
experimental value of the density of particles in liquid "^He, n — 2.17 • 10^^ cm~'^ (at no = 9%n). 
On the other hand, the density ni of supracondensate particles, calculated from Eqs. H2t)|l at T = 
at the values of the parameters indicated, is about 0.93 n, which is also in good accordance with 
experiment, taking into account that the BEG density is determined up to ±0.01 n. 

It goes without saying, the given results, which correspond to the approximation of self- 
consistent field, cannot be used directly near the A-point, where thermodynamic fluctuations 
play the key role. But the dependences, depicted in fig. 2, qualitatively correctly describe the 
temperature dependence of the density of superfluid component. 
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Fig. l.The momentum dependence of the temperature factor H82(l . obtained for different tem- 
peratures. 
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Fig. 2. The temperature dependence of the density of BEC (lower curve) and total density 
of superfluid component (upper curve), obtained according (|27|l and (|28|l for next parameters 
All = 6.14X, Dii = 2.03X, Bii = 0.00018X, A12 = 6.21K, D12 = 3.12K, B12 = 0.00225/^. The 
empirical temperature dependence of of superfluid component [Hj (circles). 
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